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Using the spectacular new high redshift supernova observations from the HST/GOODS program 
and previous supernova, CMB and galaxy clustering data, we make the most accurate measurements 
to date of the dark energy density px as a function of cosmic time, constraining it in a rather 
model-independent way, assuming a flat universe. We find that Einstein's vanilla scenario where 
px{z) is constant remains consistent with these new tight constraints, and that a Big Crunch or 
Big Rip is more than 50 gigayears away for a broader class of models allowing such cataclysmic 
events. We discuss popular pitfalls and hidden priors: parametrizing the equation-of-state wx{z) 
assumes positive dark energy density and no Big Crunch, and the popular parametrization wx{z) = 
wo + WqZ has nominally strong constraints from CMB merely because w'q > implies an unphysical 
exponential blow-up px oc e^™o^. 



The nature of dark energy has emerged as one of the 
deepest mysteries in physics. When strong evidence for 
its existence first appeared from supernova observations 
in 1998 [1, 2], the most pressing question was whether 
it was real or an observational artifact. Since then, 
the supernova evidence has both withstood the test of 
time and strengthened [3-5], and two other lines of ev- 
idence have independently led to the same conclusion: 
measurements of cosmological clustering with the cosmic 
microwave background (CMB) and large-scale structure 
(LSS) {e.g., [6, 7]) and observation of CMB/LSS corre- 
lations due to the late integrated Sachs- Wolfe effect [8]. 
Now that its current density has been accurately mea- 
sured (WMAP-HSDSS gives px{0) = (4.8 ± 1.2) x lO'^^ 
kg/m^ [7], corresponding to (9.3±2.3) x 10^^'^'* in Planck 
units and ^Ia w 0.7), the next pressing question is clearly 
whether its density px stays constant over time (like Ein- 
stein's cosmological constant) or varies. The latter is is 
predicted by most models attempting to explain dark en- 
ergy either as a dynamic substance, "quintessence" {e.g., 
[9]), or via some form of modified gravitational theory, 
perhaps related to extra dimensions or string physics 
{e.g., [10])). See [11] for reviews with more complete 
lists of references. 

The recent discovery of 16 Type la supernovae (SNe la) 
[5] with the Hubble Space Telescope during the GOODS 
ACS Treasury survey bears directly on this question. By 
discovering 6 out of the 7 highest-redshift SNe la known, 
all at z > 1.25, this search team [5] was able to pin- 
point for the first time the transition epoch from matter 
domination to dark energy domination when the cosmic 
expansion began to accelerate. It is therefore timely to 
revisit this question of if and how the dark energy den- 
sity varies with time. This is the goal of the present 
paper. Given our profound lack of understanding of dark 
energy and the profusion of theoretical models in the re- 
cent literature, we focus on measuring the function px{z) 
in as model-independent a fashion as possible, emphasiz- 
ing what we do and do not know given various assump- 
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FIG. 1: lo" constraints on the density of matter and dark en- 
ergy from SN la (Riess sample, flux-averaged with Az = 0.05), 
CMB and LSS data, all in units of the current dark energy den- 
sity. From inside out, the four nested dark energy constraints are 
for models making increasingly strong assumptions, corresponding, 
respectively, to the 4-parameter spline, the 3-parameter spline, the 
2-parameter (fcayUii) case and the 1-parameter constant w case 
(hatched). The Universe starts accelerating when the total density 
slope din p/dln(l + z) > —2, which roughly corresponds to when 
dark energy begins to dominate, i.e., to where the matter and dark 
energy bands cross. In the distant future, the Universe recoUapses 
if the dark energy density px goes negative and ends in a "Big 
Rip" if it keeps growing (din px /dln{l + z) < 0). 



tions about how px{z) is parametrized, about data sets 
used and about modeling thereof. We will see that the 
new data are powerful enough to make previous measure- 
ments of px{z) [e.g., [12-14]) tighter and more robust 
and also to extend them back firmly into the epoch of 
cosmic deceleration. 

Analysis Technique: We wish to measure the di- 
mensionless dark energy Junction, X{z) = px{z)/ pxiO), 
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the dark energy density in units of its present value. We 
do this as described in [13], fitting to SN la, CMB and 
LSS information, obtaining the results shown in Figure 
1. 

The measured distancc-redshift relations of SNe la pro- 
vide the foundation for probing the dark energy func- 
tion X{z). In a flat Universe, the dimensionless lumi- 
nosity distance dL{z)Ho/c = (1 + z)T{z), where r(z) = 
Jq dz'/E{z') is the dimensionless comoving distance and 

Eiz) = [n^^ii + zf + (1 - a,rdx{z)] (1) 

is the cosmic expansion rate relative to its present value. 
We use the "gold" set of 157 SNe la published by Riess 
et al. in [5] and analyze it using flux-averaging statistics 
[13, 16] to reduce bias due to weak gravitational lens- 
ing by intervening matter. We assume spatial flatness 
as motivated by inflation and discuss the importance of 
this and other assumptions below. We use CMB and LSS 
data to help break the degeneracy between the dark en- 
ergy function X{z) and 0„i. For the CMB, we use only 
the measurement of the CMB shift parameter [18], R = 
n^rl'^T{zcMB) = 1.716 ±0.062 from CMB (WMAP, CBI, 
ACBAR) [6, 17], where zcmb = 1089. The only large- 
scale structure information we use is the linear growth 
rate f{z2M) = 0.51 ± 0.11 measured by the 2dF galaxy 
redshift survey (2dFGRS) [3, 19], where Z2di = 0.15 is the 
eff'ective redshift of this survey and / = {dh\D / dh\a) is 
determined by solving the equation for the linear growth 
rate D, D"{t) + 2E{z)D'{t) - \^ra'^\ -V zfD = 0, where 
primes denote d/d{Hot). Note that the CMB and LSS 
measurements wc use (R and /) do not depend on the 
Hubble parameter Hq, and are quite insensitive to as- 
sumptions made about X{z). The SN la measurements 
used are also independent of Hq, since we marginalize 
them over the intrinsic SN la luminosity calibration. 

We run a Monte Carlo Markov Chain (MCMC) based 
on the MCMC engine of [20] to obtain a few million sam- 
ples of fim and X{z). The dark energy bands in Figure 
1 correspond to the central 68% of the X-values at each 
z and the matter band does the same for pm,{z)/px{0) = 

{i + zfn^/{i-n^). 

Results: Figure 1 shows our main results, the 
constraints on the dark energy function X{z) = 
Px{z)/px{0) for four difii'erent parametrizations, and il- 
lustrates that the assumptions one makes about the curve 
X(z) have an important effect on the results. The most 
common way of measuring dark energy properties in 
the literature has been to parametrize the dark energy 
function X by merely one or two free parameters, con- 
straining these by fitting to observed data. Table 1 in- 
cludes the historically most popular parametrizations, 
expressed as functions of the dimensionless cosmic scale 
factor a = (I + z)^^ . Parametrization A simply as- 
sumes that A (a) is a power law, with the single equation- 
of-state parameter w determining its logarithmic slope. 
From the identity d\n px /dlna = —3(1 + Wj), it fol- 
lows that parametrization B corresponds to the popular 



parametrization Wx{z) = wq + w'qZ [26], which has been 
widely used in the literature. It has the drawback of 
being rather unphysical for Wq > 0, with the dark en- 
ergy density px{z) blowing up as e'^"'"^ at high redshift. 
Parametrization C avoids this [21], and corresponds to 
Wx = w\+Wa{^—a), but blows up exponentially in the fu- 
ture as a — > oo for Wa > 0. In contrast, our parametriza- 
tion D remains well-behaved at all times: both early on 
and in the distant future, the dark energy approaches ei- 
ther a constant equation of state or a constant density, 
depending on the sign of {1 +Wi). 

Obviously, the more restrictive the assumptions about 
X are, the stronger the nominal constraints will be, so it 
is crucial to be clear on what these assumptions are. For 
instance, Table 1 shows that parametrizations A, B and 
C all tacitly assume that X{z) > 0, i.e., that the dark 
energy density cannot be negative, hence ruling out by 
fiat the possibility that the Universe can recoUapse in a 
Big Crunch. Note that even arbitrary function w{z) has 
this hidden assumption built in. 

To introduce as little theoretical bias as possible into 
our measurement, wc \mc parametrizations E and F 
from Table 1; these arc fairly model-independent re- 
constructions of the dark energy function X{z), assum- 
ing merely that X{z) is a sufficiently smooth function 
that it can be modeled with a cubic spline out to some 
redshift ^max, and by a constant-w power law there- 
after. We choose ^max to avoid sparse SN la data, and 
parametrize X by its values at TV equispaced spline points 
at Zma.x/N, 2zta3,x/N,...,Zma.x- X (z) IS matched smoothly 
on to (1 -I- ^)3(H-i"i) a,t z > Zmax- This specifies X{z) 
uniquely once we require X{z) and X'{z) to be every- 
where continuous and set A(0) = 1, A'(0) = X{zi)/zi. 
We have choose Zmax = 1-4, as there are only two SNe la 
at higher redshifts. Since X{z) is only very weakly con- 
strained beyond z > Zmax, we impose a prior of Wi > —2 
to avoid an unbounded parameter space. Changing the 
prior to Wi > —20 or changing the functional form of 
X{z) at 2; > Zmax (to an exponential, for example) has 
little impact on the reconstructed X(z). We also find our 
results to be rather robost to data details. Including the 
"silver" sample from [5] does not change our results qual- 
itatively, and replacing the CMB shift parameter we used 
(i? = 1.716 ± 0.062) by R= 1.710 ± 0.137 (from WMAP 
data alone [6]) broadens the 68% confidence envelope by 
less than 20%. 

Figure 1 also shows the constraints on the dark en- 
ergy function X(z) corresponding to parametrizations A 
and D from Table 1, imposing the priors Wi > —2 and 
/oo > for D.For comparison with the results of [5], 
we also studied parametrization B, with a weak prior 
Wq > —20 to avoid an unbounded parameter space. Note 
that MCMC tacitly assumes uniform prior on the param- 
eters, so if the parameter space is unbounded, the MCMC 
will drift off in the unbounded direction and never con- 
verge. Reparametrizing changes this implicit prior by 
the Jacobian of the transformation. Although we have 
imposed minimal priors to avoid unbounded parameter 
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Table 1: Parametrizations used for the dark energy function X = px{z)/px{0) in terms of the cosmic scale factor a = (1 + z) 



Parametrization 



n Parameters 



Definition 



A) Constant eq. of state w 1 w 



B) AfJine w(z) 

C) AfRne w{a) 

D) Forever regular 

E) 3-parameter spline 

F) 4-paj:ameter spline 



2 wo, w'q 

2 Wl, Wa 

2 Wi, foo 

3 Wi, X{zi), X{z2) 



X = a-3(i+'") 

X = (j-3(l+™l+«'<i)g3TOo(a-l) 

X = /oo + (1 - /cx,)a-3(i+"'*) 
Cubic spline in z for z < Z2, X = 



X{Z2) 



1 + z 

1 + 22 



J3{l+^^>i) 



for z > Z2 



4 Wi, X{zi), X{z2), X{z3) Cubic spline in z for z < zs, X = X{z3) ( ) for z > zs 



space where X{z) can be arbitrarily close to zero, but 
we have not imposed priors motivated by any theoretical 
model. For example, scalar-field models typically have 
X'(a) < 0, since fields usually roll down potentials, not 
up. In addition, many models prohibit the dark energy 
density from being negative. However, we do not wish to 
assume such priors, since "dark energy" could be a mani- 
festation of something completely different, like modified 
gravity [10]. 

As has been emphasized [22-24], SN la data are sen- 
sitive only to the smooth, overall shape of X{z). This is 
because the error bars on sharp features on a scale Az are 
proportional to (A2;)~^/^ due to the derivative involved 
in going from comoving distance r{z) to dark energy func- 
tion X [23] — reconstructing wx{z) is still harder, the 
requirement that one effectively take the second deriva- 
tive of noisy data [14] giving the error scaling as {Az)~^^^ 
[23]. Figure 1 shows that as we allow more small-scale 
freedom by parametrizing X[z) by 1, 2, 3 and 4 param- 
eters, the allowed bands become thicker. However, the 
broader bands generally encompass the narrower ones, 
showing no hint in the data that the true X{z) has funny 
features outside of the 1- and 2-parameter model fami- 
lies. Indeed, all bands are seen to be consistent with 
the simplest model of all: the zero-parameter "vanilla" 
model X{z) — 1 corresponding to Einstein's cosmological 
constant. 

In other words, faced with the fact that an analysis 
using parametrization A implies w w — 1 (we obtain 
w = -OMtoAl combining SN la, CMB and LSS), read- 
ers hoping for something more interesting than vanilla 
may correctly argue that these constraints are dominated 
by accurate measurements at lower rcdshift and may fail 
to reveal hints of an upturn in X{z) at z > 1 because 
parametrization A incorrectly assumes that (logo, logX) 
is a straight line. Our more general parametrizations 
close this loophole by allowing X{z) much greater free- 
dom, and the fact that none of them provide any hint 
yet of non-vanilla dark energy behavior therefore sub- 
stantially strengthens the case for a simple cosmological 
constant, X{z) = 1. 

What is the ultimate fate of the Universe? If for any 
of our models px eventually goes negative so that total 
density drops to zero at some time iturn, then the ex- 
pansion reverses and a Big Crunch occurs at t = 2tturn 
— this applies only if X is uniquely determined by the 
cosmic scale factor (equivalcntly z) as in Table 1, and 
not for many scalar field models [27]. The cosmic time 



t = J da/ a = J H~^dhia, and if this asymptotes to a 
finite value as a — > oo, then a cataclysmic Big Rip [15] 
occurs at this time. This is equivalent to w{z) < —1 at 
z = — 1, so parametrizations A, B and C rip ii w < — 1, 
wq — w'q < —1 and Wa > 0, respectively. 

Predictions for the future need to be taken with a 
large grain of salt, since they are obviously highly model- 
dependent. For instance, parametrizations A, B and C 
cannot crunch, whereas E and F cannot rip. Simply com- 
bining all MCMC models from all our parametrizations, 
we find that 95% of them last at least another 49 gi- 
gayears, 25% ending in a Big Crunch, 8% ending in a 
Big Rip and 67% quietly expanding forever. 

Caveats and potential pitfalls: When interpret- 
ing dark energy constraints such as those that we have 
presented, two crucial caveats must be borne in mind: 
potential SN la systematic errors and potential false as- 
sumption about other physics. We refer the reader to 
[3, 5] for thorough discussions of the former and focus on 
the latter. 

The SN la, CMB and LSS measurements we have used 
involve only X{z), and iltoi- Because of degeneracies 
between these three quantities, the inferences about X{z) 
therefore depend strongly on the assumptions about the 
two cosmological parameters Vim and Qtot- Yet it is all 
too common to constrain dark energy properties using 
prior information about flm and Otot that in turn hinges 
on assumptions about the dark energy, usually the vanilla 
assumption X{z) = 1, a pitfall emphasized by, e.g., [24]. 

We have assumed flat space, fitot = 1, as have virtually 
all recent publications measuring dark energy properties 
(usually using parametrizations A, B or C). It is well- 
known that this assumption is crucial: introducing f^tot 
as a free parameter to be marginalized over has such a 
dramatic effect on luminosity distances that essentially 
no interesting constraints can be placed on X{z) at the 
present time, not even assuming the highly restrictive 
parametrization A. We will present a detailed investiga- 
tion of dark energy independent constraints on ^}^ot from 
CMB and LSS elsewhere. 

We now turn to the issue of dark-energy independent 
constraints on flm- As emphasized by [24], assump- 
tions about flm make a crucial difference as well. As 
an example. Figure 2 shows the constaints on {wq.w'q) 
for parametrization B. The left panel illustrates that the 
constraints from SN la alone are much weaker than those 
obtained by imposing a strong prior i}„i = 0.27 ± 0.04 as 
was done in Figure 10 of [5]. Although this prior co- 
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FIG. 2: How constraints on wq and uiq depend on assumptions and 
data used. Darker shaded regions are ruled out at 95% confidence 
by SNe la alone; lighter shaded regions are ruled out when adding 
other information as indicated. 68% contours are dotted. Models 
above the dotted line end in a Big Rip. The 157 SNe la (Riess 
sample) have been flux-averaged with Az = 0.05. 

incides with the measurement of flm from WMAP and 
2dFGRS [6] , it should not be used here since it assumes 
X{z) = 1. The right panel of Figure 2 shows the ef- 
fect of including CMB information self-consistently (via 
the i?-parameter) in our constraints. We see that wq- 
values as low as —3 remain allowed, as expected given 
the above-mentioned weak Jim-constraints, and that ad- 
ditional information (in this case from LSS) is needed to 
tighten things up. This panel also illustrates the hazard 
of poor dark energy parametrizations: the seemingly im- 
pressive upper limit on Wq tells us nothing whatsoever 
about dark energy properties via SN la, but merely re- 
flects that the unphysical exponential blowup X oc e^'"o^ 



would violate the CMB constraint. 

Conclusions: In conclusion, we have reported the 
most accurate measurements to date of the dark energy 
density px as a function of time, assuming a flat universe. 
We have found that in spite of their constraining power, 
the spectacular new high-z supernova measurements of 
[5] provide no hints of departures from the vanilla model 
corresponding to Einstein's cosmological constant. This 
is good news in the sense of simplifying the rest of cos- 
mology, but dims the prospects that nature will give us 
quantitative clues about the true nature of dark energy 
by revealing non-vanilla behavior. The apparent con- 
stancy of Px {z) also makes attempts to explain away dark 
energy by blaming systematic errors appear increasingly 
contrived, further strengthening the evidence that dark 
energy is real and hence a worthy subject of study. Fu- 
ture experiments [25] can dramatically shrink the error 
bars in Figure 1, and therefore hold great promise for 
illuminating the nature of dark energy. 

Public software: A Fortran code that uses flux- 
averaging statistics to compute the likelihood of an ar- 
bitrary dark energy model (given the SN la data from 
[5]) can be found at http : / /www. nhn.ou.edu/ ~ 
wang/SNcode/. 
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